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Abstract 



We have investigated the role of the cyclotron emission associated to cosmic 
magnetic fields on the evolution of CMB spectral distortions by considering the 
contributions by spontaneous and stimulated emission and by absorption in the 
computation of the photon and energy injection rates. These cyclotron emission 
rates have been numerically compared with those of the relevant radiative processes 
operating in the cosmic plasma, bremsstrahlung and double Compton scattering, for 
realistic CMB distorted spectra at early and late epochs. For reasonable magnetic 
field strengths we find that the cyclotron emission contribution is much smaller than 
the bremsstrahlung and double Compton contributions, because of their different 
frequency locations and the high bremsstrahlung and double Compton efficiency to 
keep the long wavelength region of the CMB spectrum close to a blackbody (at elec- 
tron temperature) during the formation of the spectral distortion. Differently from 
previous analyses, we find that for a very large set of dissipation mechanisms the role 
of cyclotron emission in the evolution of CMB spectral distortions is negligible and, 
in particular, it cannot re-establish a blackbody spectrum after the generation of a 
realistic early distortion. The constraints on the energy exchanges at various cosmic 
times can be then derived, under quite general assumptions, by considering only 
Compton scattering, bremsstrahlung, and double Compton, other than, obviously, 
the considered dissipation process. Finally, upper limits to the CMB polarization 
degree induced by cyclotron emission have been estimated. 

Key words: Cosmic microwave background; Cosmic magnetic fields; Radiation 
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mechanisms; Radiative transfer; Scattering 
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1 Introduction 



The cosmic microwave background (CMB) spectrum emerges from the ther- 
malization redshift, z t herm, with a shape very close to a Planckian one, owing 
to the tight coupling between radiation and matter through Compton scatter- 
ing and photon production/absorption processes. Bremsstrahlung and double 
(or radiative) Compton were extremely efficient at early times and able to 
re-establish a blackbody (BB) spectrum from a perturbed one on timescales 
much shorter than the expansion time, t exp = a/ (da/dt), a(t) being the cosmic 
scale factor and t the time (see, e.g., Danese and De Zotti (1977)). Consider- 
ing the effect of these processes combined to that of Compton scattering, the 
value of Ztherm (— 10 6 — 10 7 ) (Burigana et al., 1991) depends on the baryon 
density parameter, Q b , and the Hubble constant, H , through the product 
n b = ^[tfo/^OKms^Mpc- 1 )] 2 . 

On the other hand, physical processes occurring at redshifts z < z t herm ma y 
lead imprints on the CMB spectrum. Therefore, the CMB spectrum car- 
ries crucial informations on physical processes occurring during early cosmic 
epochs (see, e.g., Danese and Burigana (1993) and references therein) and the 
comparison between models of CMB spectral distortions and CMB absolute 
temperature measures can constrain the physical parameters of the considered 
dissipation processes. 

In the presence of a cosmic magnetic field, another photon production/absorption 
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process, the cyclotron emission, operates in the cosmic plasma. The cyclotron 
emission (or synchrotron emission in the case of relativistic particles) could be 
polarized (Rybicki and Lightman, 1979) and its degree of polarization is an 
important indicator of the field's uniformity and structure (Widrow, 2002). 

The contribution of the cyclotron emission to the evolution of CMB spec- 
trum depends on the amplitude of magnetic field and on the electron density. 
Previous studies considered the cyclotron emission by including only the spon- 
taneous emission term (Puy and Peter, 1999) or by taking into account also 
the absorption and stimulated emission terms but assuming approximations 
for CMB distorted spectra that do not fully characterize the CMB spectral 
shapes realistically predicted in the presence of energy dissipation processes, in 
particular at long wavelengths where the cyclotron emission occurs (Afshordi, 
2002). 

In this work we derived the contribution of the cyclotron emission to the 
evolution of the CMB photon occupation number, r], as a further term in the 
Kompaneets equation (Kompaneets, 1956) by exploiting the method described 
by Afshordi (2002) to take into account the cyclotron spontaneous emission, 
absorption and stimulated emission terms and generalizing it to be able to ex- 
haustively treat various reasonable choices for the photon occupation number. 
We then apply this result to realistic assumptions for the CMB distorted spec- 
tra in order to provide robust estimates of the global photon production rate 
as a function of the relevant parameters and discuss the role of the cyclotron 
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emission in the thermalization and evolution of CMB spectral distortions. 

In Sect. 2 we briefly report on the main observational and theoretical aspects 
of cosmic magnetic fields relevant for the present work. The contribution of 
the cyclotron emission associated to cosmic magnetic fields to the evolution 
of the CMB spectrum is derived in Sect. 3, where the cyclotron frequency is 
compared with the other characteristic frequencies relevant in this context. In 
Sect. 4 we compare the production rates of photon number and energy densi- 
ties from the cyclotron emission with those from bremsstrahlung and double 
Compton for two different realistic cases of distorted spectra, Bose-Einstein 
like and Comptonization like distortions. The typical values of the cosmic 
magnetic field for which the cyclotron emission contribution to thermaliza- 
tion process could be comparable to the contribution from bremsstrahlung 
and double Compton are derived. Finally, we discuss our results and draw the 
main conclusions in Sect. 5. 

2 Magnetic fields in cosmology 

Magnetic fields pervade the universe at different scales (see, e.g., Vallee (2004), 
Gaensler et al. (2004) for recent reviews), from the scale of planet and stars 
to the scales of galaxies and clusters of galaxies and of the whole universe, 
possibly affecting the cosmogonic process (see, e.g., Subramanian and Barrow 
(1998), Rees (2000)). In this section we briefly report on the main observational 
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results on cosmic magnetic fields and on physical models for their generation 
in the early universe and their evolution (see, e.g., Kronberg (1994), Grasso 
and Rubinstein (2001), Carilli and Taylor (2002), Widrow (2002) for reviews). 

2. 1 Observations of magnetic fields 

The main observational traces of Galactic and extragalactic magnetic fields 
are the Zeeman splitting of spectral lines, the intensity and polarization of 
the synchrotron emission from free relativistic electrons, and the Faraday ro- 
tation measurements (RM) of the polarized electromagnetic radiation passing 
through a ionized medium. 

The Zeeman splitting, though direct, is typically too small to be usefully used 
for observations outside the Galaxy. 

The RM of the radiation emitted by a source at a redshift z s is given by 
(Grasso and Rubinstein, 2001, Widrow, 2002) 

RM{z s ) -8.1 x 10 5 Jj^l?l( z )(l + z )-*dl(z)^, (1) 
J cm 6 /iG ur 
o ^ 

where B\\{z) is the magnetic field strength along the line of sight and 

dl{z) ~ -^(1 + z)- l [n m (l + zf + (1 - Q m - Q A )(1 + zf + n A ]-^ 2 dz;(2) 

-no 

here Q m and are the matter and cosmological constant (or dark energy) 
density parameter, and n e is the electron number density along the considered 
line of sight. 
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The interstellar magnetic field in the Milky Way has been determined using 
several methods which allowed to obtain valuable information about the am- 
plitude and spatial structure of the field. The average field strength is found to 
be ~ 3 — 4 fj,G. Such a strength corresponds to an approximate energy equipar- 
tition between the magnetic field, the cosmic rays confined in the Galaxy, and 
the small-scale turbulent motion (Kronberg, 1994) 



The field keeps its orientation on scales of the order of few Kpc, compara- 
ble with the Galactic size, and two reversals have been observed between the 
Galactic arms, suggesting that the Galaxy field morphology may be symmet- 
rical. Magnetic fields of similar intensity have been observed in a number of 
other spiral galaxies. 

Observations on a large number of Abell clusters (Kim et al., 1991), some 
of which have a measured X-ray emission, give valuable informations about 
magnetic fields in clusters of galaxies. The magnetic field strength in the in- 
tracluster medium (ICM) is well described by the phenomenological equation 



where L is the reversal field length. Typical values of L are 10 — 100 Kpc 
corresponding to field amplitudes of 1 — 10 /iG. 

High resolution RMs towards very distant quasars have allowed to probe mag- 



B 2 



(3) 
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netic fields in the distant past. The measurements are consistent with an aver- 
age field strength of 0.4 — 4 /xG on a coherence length of ~ 15Kpc, comparable 
with a typical galaxy size. 

The RMs towards distant quasars are also used to constrain the intensity of 
magnetic field in the intergalactic medium (IGM). Assuming an aligned cosmic 
magnetic field, the RMs of distant quasars imply B IGM ~ 1CT 11 G (for a simple 
cosmological model with Vt tot = Vt m = 1, — 0, and H = 75Kms _1 Mpc _1 ). 
Unfortunately the largest reversal scale is at most ~ 1 Mpc. By adopt- 
ing this scale and applying Eq. (1) the limits on Bjgm are less stringent, 
Bigm ~ 10~ 9 G at present time. 

2.2 Generation of early magnetic fields 

Quantum field theory provides a large numbers of possible physical mecha- 
nisms that may generate magnetic fields in the early universe. 
Magnetogenesis typically requires a non-thermal equilibrium and a macro- 
scopic parity violation. This condition could have been satisfied during the 
phase transitions (PT) in the early stages of the universe. 
Some authors have shown that magnetogenesis is possible during phase tran- 
sitions in the quantum chromodynamics (QCD) era. During the quark-hadron 
phase transition a baryon-excess build up in front of the bubble wall, just 
as a consequence of the difference of the baryon masses in the quark and 
hadron phases (Cheng and Olinto, 1994). In this scenario magnetic fields are 
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generated by the peculiar motion of the electric dipoles which arises from con- 
vective transfer of the latent heat released by the expanding bubble walls. The 
field strength at the QCDPT time has been estimated by the authors to be 



The magnetic field on scales L ^> I can be estimated by performing a proper 
volume average of the fields produced by a large number of magnetic dipoles 
of size / randomly oriented in space (Hogan, 1983). Such an average gives 



Using Eq. (5) the strength of the magnetic field on the galactic length scale 
at the present time is found to be B(kpc) ~ 10 _20 G. 

Some of the ingredients which may give rise to magnetogenesis may also be 
found at the electroweak phase transitions (EWPT). 

Strong magnetic fields can be generated by a first order EWPT via dynamo 
mechanism (Baym et al., 1996). The authors estimated the average magnetic 
field strength at the present time B(R ~ 10 9 AU) ~ 1(T 17 - 1(T 20 G. 

Finally, the existence of a magnetic field at decoupling may induce a Faraday 
rotation in the CMB polarization signal. For example, for a field strength of 
~ 1CT 9 G Kosowsky and Loeb (1996) derived a rotation of ~ 1° at ~ 30 GHz 
possibly observable by future CMB polarization experiments. 



Bqcd — 10 s G on a maximal coherence length l coh ~ H ( 



-i 



QCD- 




(5) 
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2.3 Evolution of cosmic magnetic fields 



The time evolution of a magnetic field in a conducting medium is described 
by the equation (Jackson, 1978) 

f = Vx(vxB) + A V *B, (6, 

where a is the electric conductivity. Neglecting fluid velocity v, this equa- 
tion reduces to the diffusion equation which implies that an initial magnetic 
configuration will decay away in a time 

. r , AiraL 2 

Tdiff{L) = , (7) 

where L is the characteristic length scale of the spatial variation of B. In a 
cosmological framework, this means that a magnetic configuration with co- 
herence length L will survive until the present time to only if r (-^o) > ^o- In 
our convention L corresponds to the present time length scale determined by 
the Hubble law 



L^L(t,) a M, (8) 
a{U) 



where a(t) is the cosmic scale factor and L(ti) is the length scale at the time of 
the formation of the magnetic configuration. As shown by Eq. (7) the relevant 
quantity controlling r di ff is the electric conductivity of the medium. This 
quantity changes in time depending on the varying population of the available 
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charge carriers and on their kinetic energy. Assuming that all charge carriers 
are non-relativistic, the estimate of magnetic diffusion length is simple. For 
simplicity, we consider only one charge carrier type, the electrons, with charge 
e, number density n e , mass m e , and velocity v. Comparing the Ohm law 
J = ctE with the current density definition J = n e ev, and using the expression 
v ~ eEAr/m e for the mean drift velocity in the presence of the electric field E, 
At being the average time between collisions of the considered charge carrier, 
for the electron conductivity we have: 



n e e 2 Ar 

a = . 9 



For the evolution of the electron number density n e we can use here 2 the 
usual formula 

n e ~ (7/8)n 6 ~ 2.45 x 1(T 6 Cl b (l + zf cm" 3 , (10) 



where a primordial helium abundance of ~ 25% by mass has been assumed 
here for numerical estimates (in this hypothesis nn = 0.75rib, and n# e = 
(l/16)n b , and n b = 2.8 x 10 -6 0^(1 + z) 3 cm -3 ). Since electron resistivity 
is dominated by Thomson scattering off cosmic background photons then 
At ~ l/n^a T c, where ar = 87r/3(e 2 /mc 2 ) 2 is the Thomson cross section, 

2 In this work, we are focussing on the pre-recombination era when matter were 
highly ionized. 
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and therefore Eq. (9) gives 



a = 



m e aTn 7 c 



(11) 



where n 7 ~ n P ~ 4 x 10 2 (T /2.7K) 3 (1 + z) 3 cm 3 is the photon number 
density in the blackbody limit, aT being the present CMB energy density 
(T ~ (2.725 ± 0.002)K; Mather et al. (1999)). The high conductivity of the 
cosmic medium has a relevant consequence for the evolution of magnetic fields. 
The magnetic flux through any loop moving with fluid is a conserved quantity 
in the limit a — > oo. In fact, the diffusion equation (6) after few vector algebra 
operations implies 



^#2 = f ^_Vx(vxB)rfS = -^- f VxVxB.rfS, (12) 
dt J dt v 1 Ana J v ' 

S(t) S(t) 



where S is any surface delimited by the loop. On scale where diffusion can be 
neglected the field is said to be frozen-in, in the sense that lines of force move 
together with the fluid. Assuming that the universe expands isotropically, 
magnetic flux conservation implies 



2 



B(t) = B(t ) (^J =B (1 + .) 2 , (13) 



where B is the present time magnetic field. 
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3 Effect of the cyclotron emission associated to cosmic magnetic 
fields on the evolution of CMB spectral distortions 

The evolution of the CMB photon occupation number, rj(u,t), at redshifts 
z £ 10 6 — 10 7 , is well described by the Kompaneets equation (Kompaneets, 
1956). We can write formally this equation as 



where A, take into account processes that do not change the photon number 
and Tj take into account photon production/absorption processes. 
During the cosmic epochs of interest here, before recombination, the relevant 
processes are Compton scattering, bremsstrahlung (BR) and double (or ra- 
diative) Compton (DC). We want to study the role of magnetic fields on the 
evolution of the CMB photon occupation number. The effect of a magnetic 
field in a ionized plasma is to speed up all the present charged particles. We 
consider here only the electrons. 

We shall show in this section how to calculate the contribution (drj/dt)cE 
in Eq. (14) due to the cyclotron emission (CE) of electrons accelerated by a 
cosmic magnetic field in the primeval plasma. Then we compare this contri- 
bution with that of bremsstrahlung and double Compton in the case of early 
and late CMB spectral distortions by using appropriate analytical formulas 
for the description of the photon occupation number 77. 




(14) 
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3.1 Cyclotron emission 



Following the approach presented by Afshordi (2002) the rate of energy loss 
via cyclotron emission by non-relativistic electrons moving in a magnetic field 
B can be obtained classically (Jackson, 1978) 

d£ 2e 2 u 2 (v ± 2 )n e 2 e 4 B 2 (v ± 2 )n e 



dt 3 c 3 3 m 2 & 



(15) 



here uj c = 2nu c = eB/m e c is the cyclotron frequency, E c = hu c = h\v c is the 
energy of the emitted photon, v± is the component of the electron velocity 
normal to the magnetic field direction. In the non relativistic limit, almost 
all the emitted photons have the frequency uj c and thus the rate of photon 
production per unit volume, ip, can be obtained using Eq. (15) 



_ dS/dt _ 2 n e e 3 B(v ± 2 } 
hu r 3 hm P c 4 



By assuming a Maxwellian distribution for the electrons we have 



(v± 2 ) = ltf) = ^, (17) 



where T e is the temperature of the electron gas. Eq. (16) then becomes 
An e e 3 Bk B T e 

^=3 Km 2 c 4 ■ (18) 



The presence of photons in the environment enhances the photon production 
through stimulated emission. Also photons can be absorbed by the rotating 
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electrons. These processes can be expressed via (Spitzer, 1998) 




) = J2 A (^ + v(E'))Ve(E + E c )-Br ] (E')r ]e (E) 5{E' - E c ), (19) 



/ CE [{E} 



where r\ e is the electron distribution, A e B are the Einstein coefficients, E' is 
the photon energy, E + E c is the electrons energy and S(E' — E c ) is the Dirac 
5 function. The sum is over the energy states of the electrons (Landau levels). 
The first two terms in the second member of Eq. (19) describe the sponta- 
neous and stimulated emission, the third term is the contribution due to the 
absorption of photons: a photon of energy E' = E c can be emitted by an 
electron which undergoes a decrease of energy from E + E c to E, or may be 
absorbed by an electron of energy E (which clearly raises its energy to E + E c ). 
The coefficient B can be obtained considering that for a Planck distribution 
drj/dt = 0; therefore 



By assuming a Maxwellian distribution for the electrons, in the non relativistic 
limit we can write 



£ eW» T lAn,(E + E c ) = £ B Ve (E) . 



(20) 





from which one obtains 




From these equations together with Eq. (20) we have 
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B = A e (E ~ Ec)/kBTe . 
Substituting this expression in Eq. (19) is straightforward to obtain 

= An e e- E ^ kBT ^[l + rj(E') - e E ' l kBT ^{E')]8{E' - E c ) . (21) 

\ ^ / CE 

We can integrate Eq. (21) over the phase space to have the total photon 
injection rate 

here the factor 2 takes into account the possible polarization of the photon 
and the coefficient A in Eq. (21) is the photon production rate for zero photon 
occupation number, i.e. the same as ip in Eq. (18); inserting the expression for 
drj/dt given by Eq. (21) in Eq. (22) one finds 

A ^^eck B T e eEc/kBTe _ (23) 

By introducing the dimensionless frequency x e = E' /ksTe, from Eqs. (23) and 
(21) we have 

= K CE (z) [l-r, (x e ) (e x * - 1)] 5(x e - x e , CE ) , (24) 



where 



4tt^ 6 c * 
K CE (z) = -^—n e = 4.64xlO- A n b -B- 1 (z)(l + z) 3 s" 1 . (25) 
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The magnetic field will be assumed to scale as in Eq. (13). 
Now we can write the complete Kompaneets equation including cylotron emis- 
sion 



di] 1119 

dt (ft tc x 2 dx 



x 



+ 



K 



BR- 



9BR 



x- 



Kdc^y + K CE5{x e - X e>C E) 



(26) 

[1-^.-1)] , 



where the coefficients K(z) and the Gaunt factors, qbr and goc, for bremsstrahlung 
(Karsaz and Latter, 1961, Rybicki and Lightman, 1979) and double Compton 
(Gould, 1984) are given in Burigana et al. (1991) and Burigana et al. (1995), 
tc = mc 2 /[kT e (n e aTc)] is the timescale for the achievement of kinetic equilib- 
rium between radiation and matter, <p = T e /T r where T r = T (l + z) is the 
CMB temperature, and x = (px e . 

At this point a little consideration on the cyclotron term is necessary. In a real- 
istic framework cyclotron emission should not be a line emission, as considered 
here, but it should be a continuum emission peaked around the characteristic 
frequency vqe- In addition, fluctuations in the plasma and in the magnetic 
field imply analogous fluctuations in the cyclotron emission frequency. There- 
fore a better approximation should involve a replacement of the 5 function 
approximation with a proper (physically motivated) smooth function in all 
numerical estimates. On the other hand, as it will be evident from the results 
presented in the next sections, the simple 5 function approximation does not 
significantly limit the validity of all our main conclusions that rely on the 
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typical magnitude of the cyclotron emission frequency as function of redshift 
and not on its detailed value. 



3.2 Numerical comparison between characteristic frequencies 



In this section we numerically compare the characteristic frequencies appear- 
ing in the "classical" (i.e. in the presence of C, DC, and BR) solutions of 
the Kompaneets equation with that of cyclotron emission. We can write the 
dimensionless frequency of cyclotron emission as follow: 



X e ,CE 



hVc ~ 4.97 x 10- 5 • 0- 1 ( ) B Q (1 + z), 



kT P 



2.7K 



(27) 



where B is the present value of the magnetic field. Let us define the frequency 
the solution of the equation t a i> s = tc (Sunyaev and Zeldovich, 1970, 
Illarionov and Sunyaev, 1974), where t abs is defined by (Burigana et al., 1991, 
Burigana et al., 1995) 



l/t abs = E BR (x e ,z)e Xe +E DC (x e ,z) (e Xe -l) 



K B R(z)g B R(x e )e Xe K DC {z)g D c{x e ) 



+ 



xvl 



(e- - 1) 



(28) 



Let us define also the frequency x e ^ s according to the relation (Danese and 
De Zotti, 1980): 



l+z h 

y a bs(xe,t,t h ) = r(t) - r{t h ) = J 

1 + 2 



t exp d(l + z) 

tabs 1 + Z 



(29) 
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when i/abs is evaluated at the current time to (z 
(redshift) at which a given distortion occurred. 



= 0); here th (zh) is the time 




7 6 5 4 3 

Log( 1 +z) 



Fig. 1. Comparison between the characteristic frequencies x e)C (thick solid line), 
x e ,abs (thin solid line) and x ej cE, for different values of the magnetic field Bq: 10 _9 G 
(thick dashed line), 10~ 8 G (thick dott-dashed line), 10~ 7 G (thick three dots and 
dash line). In the three panels different values of the fractional energy Ae/ej, injected 
in the CMB radiation field are assumed. We consider the case of early (Bose-Einstein 
like) distorted spectra with Ae/ei (~ no/lA for fiQ <C 1) values corresponding to 
the following values of the chemical potential (see also Sect. 4.2): I0~ 4 (panel a), 
10~ 2 (panel 6), 1 (panel c). The cosmological parameters relevant in this context 
have been assumed in agreement with WMAP (see Tab. 3 of Bennett et al. (2003)); 
nevertheless the dependence on the detailed choice of the assumed cosmological 
parameters is here not particularly critical. 
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The frequency x eta b s is the maximum frequency at which the considered photon 
emitting processes could re-estabilish a thermal (blackbody) spectrum after 
the distortion. On the other hand, for x e £ x eyC the timescale tc is shorter 
than t a b s and the high Compton scattering efficiency prevents the formation 
of a blackbody spectrum and, as a particular case, for z > Z\ leads to a 
Bose-Einstein like spectrum with a frequency dependent chemical potential 
(possibly evolving in time) tending to vanish at x e <C x e>c and to a value 
independent on the frequency at x e ^> x e>c , because of the different relative 
relevance of Compton scattering and photon emission processes. Here z\ is the 
redshift for which the dimensionless time variable y e (z) defined by 



assumes the value y(zi) ~ 5 (Zeldovich and Sunyaev, 1969, Illarionov and 
Sunyaev, 1974, Chan and Jones, 1975, Burigana et al., 1991) and the kinetic 
equilibrium between matter and radiation can be reached. 

In Fig. 1 we show the comparison between the dimensionless frequency x £! ce(z) 
and the other characteristic frequencies x eja b s (z) and x eiC (z) for different val- 
ues of the energy injection rate Ae/ei and of the magnetic field B , in the 
cosmic epochs of interest here. As evident, x €j ce <C ^e,c and x £j ce x eA b s at 
any time. This implies that, in the presence of a dissipation mechanism op- 
erating to perturb a blackbody spectrum through a non vanishing Compton 
scattering term associated to a disequilibrium between matter and radiation 




(30) 
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temperatures, the cyclotron emission process term occurs at frequencies where 
bremsstrahlung and double Compton are very efficient in keeping the radia- 
tion spectrum close to a blackbody at matter temperature. In the next section 
we will quantitatively investigate this point. 

4 Results 

4-1 Comparison between photon and energy injection rate 

In this section we compare the contributions given by the different photon 
production terms in the Kompaneets equation for various cosmic epochs. We 
will compute the photon number production rate 




(31) 



and the photon energy production rate 




(32) 



Using Eq. (24) for the CE term we find 




-jf) Kce(z) f[l-ri(x e )(e°>-im 
{k B T e f K CE (z) J [I- ri(x e )(e** - l)]<f(s 



x e — x e ^cE)x\dx e 




(33) 
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Using the expressions given by Burigana et al. (1991) and Burigana et al. 
(1995), for the BR and DC term we have 



[fj B = *»M^ /^[l - - 1)] e-d*. 

( dr\ R<rrk A T A r 

Ut) =Kbr{z) ~^^ j9BR(x e )[l-ri(x e )(e°>-l)]e-*dx e , 

\ / BR 

(34) 

(1)^^)5^/ ^-^-1)1 .^ iXt 

(§) =Knc(zf-^ J [1 -v(x e )(e°> - 1)] dx e , 

\ / DC 

where the double Compton Gaunt factor (in the limit of a BE-like spectrum; 
Burigana et al. (1995)) has been explicitely included. 

4-2 Early distortions 



We compare the contributions by the three radiative processes considered 
here for an energy dissipation occurring at high redshifts by assuming for r\ 
(i.e. in Eqs. (33) and (34)) a pure Bose-Einstein (BE) formula or a BE-like 
spectrum with a frequency dependent chemical potential (Eqs. (35) and (36)), 
as expected under realistic conditions for the active phase of a dissipation 
process and the subsequent CMB spectrum evolution (Sunyaev and Zeldovich, 
1970, Illarionov and Sunyaev, 1974): 



Vbe = +fi( * e) _ I , (35) 



fi{x e ) = fi e~ x ^ . (36) 
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Neglecting the subsequent decreasing of the chemical potential due to the pho- 
ton emission/absorption processes, the amount of fractional energy dissipated 
in the plasma for a process with a negligible photon production is related to 
fi by a well known relation (Sunyaev and Zeldovich, 1970, Danese and De 
Zotti, 1977) that in the limit of small distortions is simply expressed by 

fio ~ lAAe/ei , (37) 

where q is the radiation energy density before the beginning of the process. 

The result of our computation is reported in Figs. 2 and 3 for a representative 
choice of the cosmological parameters and some values of the cosmic magnetic 
field 3 . 

Considering the pure BE approximation, exploited for example by Afshordi 
(2002), and neglecting the consideration in footnote 3, the CE contributions 
appear to be much larger than the contributions by BR and DC even for 
B & 10~ 8 , i.e. for realistic values of magnetic fields (see panels a\, a2, h\ and 62 
in Figs. 2 and 3). Since this result is obtained under the assumption that the 
chemical potential \x is constant with respect to the frequency it neglects at 

3 Note that while all the integrals in Eq. (34) converge and do not depend on the low 
frequency cut-off, x;, adopted in the numerical integration (provided that xi <C x etC ), 
in the case of a photon distribution function described by the pure BE formula the 
photon number production rates do not converge because of the low frequency l/x e 
dependence. The results reported in panels a\ and 02 of Figs. 2 and 3 for the DC 
and BR photon number production rates, obtained by using logx; ~ —4.3, are only 
indicative. Of course, this does not affect the conclusion of our work that is based on 
the proper BE- like spectrum which frequency dependent chemical potential removes 
the above divergency (see panel a of Figs. 2 and 3). 
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Fig. 2. Photon number (panels a and a\) and associated energy injection (panels b 
and b\) fractional (i.e. divided respectively by the photon number and energy den- 
sity of a blackbody at temperature Tq(1 + z)) rates at the redshift z from DC (dots), 
BR (long dashes) and CE (dashes: B = 2 x 1(T 6 G; dots-dashes: B = 2x 1(T 5 G; 
three dots-dashes: Bq = 2 x 10~ 4 G) in the presence of an early distortion with 
Hq = 1.4 x 10~ 2 occuring exactly at the redshift z (properly speaking this result 
holds for z £ z\ while it is only indicative for BE like distorted spectra at z < z\ 
resulting from processing occurred at z ^ z±). The panels a and b refer to compu- 
tations including the proper frequency dependence of the chemical potential [i.e. a 
BE-like spectrum, fj, = /i(x e )], neglected for comparison in panels a\ and b\ (i.e. 
fi = Ho)- Panels 02 and 62 are identical to panels a\ and b\ but for Bq = 1CT 9 G 
(dashes), Bq = 10 -8 G (dots-dashes), Bq = 10 -7 G (three dots-dashes) and show 
that neglecting the frequency dependence of the chemical potential would imply the 
(wrong) conclusion of a dominance of the CE contribution for reasonable values of 
the magnetic field. Note that the DC and BR rates in panels a\ and 02 are only 
indicative (see footnote 3). The values used for the cosmological parameters are the 
same as in Fig. 1. 
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Fig. 3. The same as in Fig. 2 but for = 1-4. 



all the efficiency of bremsstrahlung and double Compton at low frequencies in 
the formation of the distorted spectrum during the energy dissipation phase. 



In reality, the result is very different when we adopt a better approximation 
of the distorted spectrum, represented by the BE-like spectrum as discussed 
above. The differences found for CE term by assuming a constant or a fre- 
quency dependent chemical potential are much larger than those found for the 
DC and BR terms (see Figs. 2 and 3) because of the different frequency loca- 
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tion of their contributions (see Fig. 1). Replacing the expression for 77 given 
in Eqs. (35) and (36) in Eqs. (33) one obtains in fact 



dn 
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Since for z > z± it is £ ei aE <^ x e ,c (see Fig. 1), from Eq. (36) it is straight- 
forward to obtain ^(x^ce) <^ 1 (f° r B & 10 -6 ), and the above rates become 
then very small or negligible (Zizzo, 2003). We find in fact that the CE term is 
comparable to the BR and DC term only for values of B larger than ~ 10~ 6 G 
(see Fig. 4 and the following discussion), the exact value depending on Ae/ei 
and on the cosmological parameters. This value of B is much larger than the 
values obtained from observations in IGM (see Sect. 2.1) and of the predictions 
of the scenarios for the genesis of cosmic magnetic fields (see Sect. 2.2). 

Differently from previous analyses, we conclude that the CE contribution to 
the thermalization process of CMB spectrum after an early heating is negligi- 
ble. 

For sake of generality, we have computed the value of the cosmic magnetic field 
for which the CE photon production rate is equal to (or is 1/10 or 1/100 of) 
the combined photon production rate by BR and DC. We report in Fig. 4 the 
results found at the two representative redshifts, z = 10 7 and z = 3 x 10 6 , quite 
close to the thermalization redshift, since the relative contribution by the CE 
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term to the global photon production is more relevant at higher redshifts (see 
Figs. 2 and 3). As evident from Fig. 4, for Ae/e^ ^ some x 1CT 2 the result is 
very weakly dependent on Ae/e« while it significantly depends on the assumed 
ratio between the CE rate and the BR and DC combined rate. Exploiting 
the results at Ae/e; <, few x 10~ 2 , we have then searched for a fit of B as a 
function of the ratio of these photon production rates in the form: 



where the overdots denote the time derivatives. The best-fit values of the 
parameters are a ~ 0.178 and (3 ~ 1.924 for z h = 10 7 and a ~ 0.252 and 
P ~ 3.064 for z h = 3x 10 6 . 

4-3 Late distortions 

In the case of small, late distortions C 1), a quite good approximation of 
the high frequency region of the CMB spectrum (observationally, at A £ 1 cm), 
where Compton scattering dominates, is provided by the well known solution 
(Zeldovich and Sunyaev, 1969) 




(38) 



T)(x,T)~r) (x,T) —T]i + U 



x/(j)jexp(x/(f)i) ( x/4>i 



-4 



(39) 



[exp(x/(f>i) — l] 2 \tanh(x/20j) 



where rji is the initial distribution function. 



Here u is the (evolving) Comptonization parameter 
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Fig. 4. Magnetic field Bq for which the ratio between the photon number production 
rates by cyclotron emission and bremsstrahlung plus double Compton assumes the 
values indicated close to each curve (1, 0.1, 0.01) as a function of the of fractional 
injected energy. The values used for the cosmological parameters are the same as in 
Fig. 1 (see also the text). 



u(t) = u(z) 



) — fa dt 

U 1 + 2 



m e c L 



1 + z' 



By integrating u over the whole relevant energy dissipation phase we have 
the "usual" Comptonization parameter related to the whole fractional energy 
exchange by the well known expression (Zeldovich et al., 1972) 



u ~ (l/4)Ae/ei 



(41) 
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In the case of an initial blackbody spectrum rji is given by 



1 



(42) 



Vbb 



where 



^ = (l + A£/ £i )- 1/4 -l 



— u 



(43) 



is the ratio between electron and radiation temperature before the beginning 
of the process (only small modifications of above formulas are needed for initial 
BE like spectra; see Burigana et al. (1995)). 

For z < mm(z p ,z), with z = min^i, Zh), z p being the redshift 4 z p = 
2.14 x 10 4 (T /2.7^) 1 / 2 ^ 1/2 (k/1.68) 1 / 4 (Burigana et al, 1991, Burigana et 
al., 1995), the CMB spectrum at very low frequencies (observationally, at 
A £ some dm) is mainly determined by the BR term (the DC efficiency sig- 
nificantly decreases at decreasing redshifts) and it can be a quite well ap- 
proximated by the expression (Danese and De Zotti, 1980, Burigana et al., 
1991, Burigana et al., 1995) 

r)(x, t) ~ r] L (x, t) 



4 k (~ 1.68 for 3 species of massless neutrinos) takes into account the contribution of 
relativistic neutrinos to the dynamics of the universe. Strictly speaking the present 
ratio of neutrino to photon energy densities, and hence the value of k, is itself 
a function of the amount of energy dissipated. The effect, however, is never very 
important and is negligible for very small distortions. 




(44) 
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where a constant value of </>(t), related to the value of the fraction of injected 
energy, Ae/e,, has been assumed for simplicity, in order to make the treatment 
fully analytical. 

At intermediate frequencies both BR and C are important and the spectrum 
is well described by the expression (Burigana et al., 1995) 

rj(x, t) ~ r] L (x, t) + r] H (x, r) - 77, . 

In Fig. 5 we compare the rates of CE with those of BR and DC using for r\ the 
expression in Eq. (39), which, as already mentioned, is a good approximation 
for the high frequency range and describes the spectrum when only Compton 
scattering operates 5 . 

However, in the frequency interval where CE is located (i.e. i e « 1) the 
spectrum is well approximated by Eq. (44). Assuming a constant value of 
4>(t) implies that we can easily perform the integral in Eq. (44) to have 



/x r / m , 1 exp[-(r-r h )] 
V (x, r) = Vi exp[-(r - r h )\ + _ _ - 

j g Vabs 

= rj e -y ab s _\ (45) 

'* e x/<t>f — I e x/<t>i _ 1 ' V 7 



where 0/ = 0(r) and Eq. (29) has been used in the last equality. 

From Fig. 1 it is clear that for z < Z\ we have x^ce{z) x e ,abs(z), then 
Uabs{xe,CE) 3> 1 (see Figs. 6 and 7) and consequently e~ Vabs — > 0; Eq. (45) 



5 Note the analogy with the result obtained in the limit of a pure BE spectrum in 
the case of early dissipation processes. 
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Fig. 5. Photon production rates (panels a and a±) and energy production rates 
(panels b and b\) for the three processes taken into account. The expression adopted 
for 7] is given by the Eq. (39). We have assumed here: Ae/ej = 10 -4 ; Bq = 1CP 9 
(short dashed line), Bo = 10~ 8 (dot dashed line), and Bq = 1CT 7 G (three dot 
dashed line); = 10 5 (panels a and b) and z/, = 3x 10 3 (panels a\ and b\). Thick 
solid line: double Compton scattering rates; dots: bremsstrahlung rates. The values 
used for the cosmological parameters are the same as in Fig. 1. 

reduces then to 



*/<t>f - 1 



Replacing this expression for rj in Eq. (24) gives 
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Fig. 6. Optical depth y a b s at the frequency x e> cE(zh) (reported close to each curve) 
evaluated by computing the integral in Eq. 29 but within 1 + z and 1 + Zh, as 
function of z. We assume Zh = 10 5 in panel a and Zh = 3 x 10 3 in panel b. Again 
B = 1CT 9 , B = 1CT 8 , and B = 1CT 7 G (from the top to the bottom). Note that 
keeping x e< cE at the value it has at the largest redshift in the interval [z, Zh] implies 
an overestimation of x et cE and then an underestimation of the photon production 
by bremsstrahlung and double Compton. The values of y a bs reported in the figure 
have then to be considered as lower limits. The values used for the cosmological 
parameters are the same as in Fig. 1. 



0. 



This result, that can be generalized to thermal histories with variable 0(r), 
is a consequence of the fact that for z < z p photons of frequency x e < x eta b s 
are absorbed before Compton scattering can efficiently shift them to higher 
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Fig. 7. Optical depth y a fe s at the frequency x e> cE(z) evaluated by computing the 
integral in Eq. (29) but within a small interval Az (~ z/100) about z, as function 
of z, for the same values of Bq as in Fig. 6. Analogously to Fig. 6, we report in 
two distinct panels the case at z between 10 5 and 10 3 (panel a) and the case at z 
between 3 x 10 3 and 10 3 (panel b; in this context, it is in practice a zoom up of 
panel a). The values used for the cosmological parameters are the same as in Fig. 1. 

frequencies during the active phase of the dissipation process, as illustrated 
by Fig. 7: even integrated over a small interval in redshift, taken sufficiently 
small to assure that 0(r) could be considered as a constant in practice within 
the interval (as possible for reasonably smooth thermal histories with variable 
4>{t))i the (bremsstrahlung dominated at these redshifts) optical depth y abs at 
the (very low) cyclotron frequency is so large to keep the spectrum very close 
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to a blackbody at electron temperature at each time. 
5 Discussion and conclusion 

In this paper we have investigated the role of the cyclotron emission associated 
to cosmic magnetic fields on the evolution of CMB spectral distortions. We 
properly included the contributions by spontaneous and stimulated emission 
and by absorption. We have derived the photon and energy injection rates due 
to cyclotron emission exploiting and generalizing the approach by Afshordi 
(2002). The photon and energy production rates by cyclotron emission have 
been numerically compared with those of the relevant radiative processes op- 
erating in the cosmic plasma, bremsstrahlung and double Compton scattering. 
We have evaluated the role of cyclotron emission by adopting realistic approx- 
imations for the CMB distorted spectra. In the case of early distorted spectra 
we adopted a Bose-Einstein distribution function with a frequency dependent 
chemical potential, as realistically expected by an early dissipation process. For 
realistic values of -Bo, consistent with observations, we find that the cyclotron 
emission term is negligible with respect to that from bremsstrahlung and dou- 
ble Compton, because of the different frequency location of their contributions 
and the very high efficiency of bremsstrahlung and double Compton to keep 
the CMB spectrum very close to a blackbody shape (at electron temperature) 
during the spectral distortion formation. 

We have also considered the role of cyclotron emission for a dissipation pro- 
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cess possibly occurred at low redshift (zh ~ Zi) by exploiting the available 
analytical description of late distorted spectra. Again the very low frequency 
location of the cyclotron emission implies that it could be relevant only where 
bremsstrahlung and double Compton very efficiently keep the blackbody shape 
of the CMB spectrum during the active phase of the dissipation process and 
the subsequent evolution. Again, the cyclotron contribution is found to be 
negligible. 

Of course, our approach does not apply to possible specific mechanisms of 
generation of spectral distortions in which the energy is injected in the cos- 
mic radiation field through a relevant modification of the photon occupation 
number at very low frequencies, close to the cyclotron emission frequency. A 
specific (typically, highly model dependent) treatment is needed in this case. 

Differently from previous analyses, we then conclude that for a very large set 
of dissipation mechanisms, as in the case of dissipation processes mediated 
by energy exchanges between matter and radiation or associated to photon 
injections at frequencies significantly different from the cyclotron emission 
frequency, the role of cyclotron emission in the evolution of CMB spectral 
distortions is negligible for any realistic value of the cosmic magnetic field. In 
particular, it cannot re-establish a blackbody spectrum after the generation of 
a realistic early distortion and it cannot produce a significant spectral distor- 
tion (and therefore no significant limits on the cosmic magnetic field strength 
can be set by constraints on CMB spectral distortions when interpreted as 
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produced by cyclotron emission). 

Consequently, the constraints on the energy dissipations at various cosmic 
times set by the currently available data (Salvaterra and Burigana, 2002) and 
expected by future experiments (Kogut, 1996, Fixsen and Mather, 2002, Buri- 
gana and Salvaterra, 2003, Burigana et al., 2004) can be derived, under quite 
general assumptions, by considering only Compton scattering, bremsstrahlung, 
and double Compton, other than, obviously, the considered dissipation pro- 
cess (es). 

As a particular application of CMB spectrum analyses to cosmic magnetic 
fields, Jedamzik et al. (2000) estimated a (weak) limit on the magnetic field 
strength B & 10~ 8 G on a coherence length of ~ 400 pc by comparing the 
constraints on CMB Bose-Einstein like distortions (attenuated by taking into 
account the dominant double Compton process 6 ) with the distortion level 
predicted by the damping of magnetic field normal oscillation modes before 
recombination, which associated dissipation rate is proportional to the square 
of the magnetic field strength. 

An interesting aspect is the evaluation of the polarization degree imprinted 
by the cyclotron emission on the CMB. While the (very small) effect on CMB 
spectral distortions shows a strong dependence only on the magnetic field 
strength, the polarization imprint depends in principle also on the unifor- 

6 Note that 10 -8 G is a value much smaller than that for which the role of cyclotron 
emission is comparable to that of double Compton scattering (see Fig. 4). 
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mity and structure of the magnetic field and, from the observational point of 
view, on the experiment angular resolution because of the obvious decreas- 
ing of the global polarization signal when it is measured within a solid an- 
gle larger than the field coherence angular scale. A detailed study, out of 
the scope of the present work, is therefore model dependent. On the other 
hand, an approximate upper limit to the polarization degree imprinted by 
the cyclotron emission on the CMB can be derived by exploiting Eq. (38) 
at the representative redshift z~3x 10 6 , i.e at an epoch significantly later 
than the thermalization redshift of an arbitrarily large spectral distortion but 
before the strong decreasing of the cyclotron photon production efficiency 
(see Fig. 1 and panels a) of Figs. 2 and 3). Assuming that after an en- 
ergy dissipation producing a distorted spectrum with a chemical potential 7 
/i ~ 1.4Ae/ej the photon production processes are able to fully re-establish 
a Planckian spectrum, the additional photon fractional number is An/n r ~ 
[n P (T r ) -n P (T r )(f) 3 BE (ii )(p(ii )]/np(T r ); since (/>beM - (1 - l.H/i )~ 1/4 and 
(fd^o) ~ 1 — 1.368^0 we obtain An/n r ~ 0.779Ae/ej. From Eq. (38) we have 
that only a fraction r ~ 6.94 x 10~ 13 x io 3 - 97 ( B o/ 1(r6G ) of the above additional 
photons is produced by cyclotron emission. Assuming an intrinsic polariza- 
tion degree f % v for the cyclotron emission from each coherence cell, a beam 
depolarization factor fp d due to the averaging of polarization signals with 
different orientations within the experiment beam size, a Faraday depolar- 

7 Only small distortions are compatible with current data at z ~ 3 x 10 6 (Salvaterra 
and Burigana, 2002). 
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ization factor /jf , and a polarization decreasing factor f* c due to the ran- 
domization of the polarization orientations when (because of Compton scat- 
tering) cyclotron emission photons migrate from frequencies ~ u c to higher 
frequencies where CMB experiments are carried out, we derive an upper limit 
of d ~ 5.41 x 10- 13 (^/ p fed /™/ p sc )(Ae/ ei ) x k^Bo/io^g) for the contribu . 
tion of cyclotron emission to the CMB polarization degree. Even assuming 
fpfp d fp d fp C ~ 1) A e / e « ~ 0.1, and B ~ 1/zG we find an (extremely generous) 
upper limit of d ~ 5 x 1CT 10 , a value much smaller than the intrinsic polariza- 
tion degree of CMB anisotropies and of the polarization anisotropies "directly" 
induced by cosmic magnetic fields (see, e.g., Barrow et al. (1997), Subrama- 
nian et al. (2003), and references therein) and well below any observational 
chance. 
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